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We investigate the occurrence of singlet bound states in the Raman spectrum of dimerized spin 
1/2 chains by Exact Diagonalization and Density Matrix Renormalization Group techniques. We 
predict that several bulk 5-peaks could be observed in pure systems. Furthermore, we show that 
new low energy lines arise from non-magnetic impurity doping. These features are interpreted in 
terms of soliton-antisoliton and soliton-impurity bound states respectively. Energies and spectral 
weights associated with these bound states are sensitive to lattice relaxation effects. Our results 
are discussed in the context of the inorganic spin-Peierls compound CuGeOs and quantitatively 
compared to recent Raman experiments. 



I. INTRODUCTION 

One-dimensional quantum spin chains have drawn 
great attention since the discovery of a spiBrPeierls tran- 
sition in the inorganic CuGeOs compoundEJ. The transi- 
tion is experimentally inferred from an isotropic drop of 
the magnetic susceptibility below a critical temperature 
signaling the opening of a singlet-triplet spin gap A 01 
together with a lattice dimerization. 

The occurrence of well-defined magnetic bound states 
below the excitation continuum has been intensively in- 
vestigated both theoretically and experimentally. A 
triplet bound state with A 01 ~ 16.8 cm^ 1 has been 
observed by inelastic neutron scattering experimentsEm. 
Raman scattering experiments have identified a sharp 
singlet resonance-at an energy A 00 ~ 29.7 cm -1 above 
the ground stateoQ. Recently, Raman experiments on 
Zn doped CuGeOs have detected a new singlet bound 
state at A™ 13 ~ 15.0 cm -1 which is close to the value 
of the spin gap A 01 (Rcf. ||). Although it is commonly 
believed that these features signal the presence of bound 
states, there is still some controversy about their inter- 
pretation. It has been proposed that the singlet bound 
state is a magnon-magnon bound stated. An alternative 
description of both the singlet and triplet bound states, 
in terms, pf soliton-antisoliton (ss) bound states, is also 
possible&tj. However, the question remains whether such 
ss bound states are observable in Raman spectroscopy. 
R is the purpose of the present paper to show that this 
is indeed the case: states which are clearly identifiable 
as ss or impurity-soliton bound states yield a significant 
Raman spectral weight at frequencies roughly coincident 
with the experimental results. Hence, a complete under- 
standing of the experimental results is possible within the 
soliton picture both in pure and doped materials. 

Quasi one-dimensional spin-Peierls compounds are 



widely described ip the literature as frustrated dimerized 
Heisenberg chainsEJO: 



H = ^[J(l + Si-lf^i.Si+r + hSi.S., 



i+2\ 



(1) 



where the explicit dimerization S arises from some cou- 
pling to the three-dimensional lattice treated at the 
mean-field level. An extension of this model allowing for 
(slow) variation of the order parameter 6 in space will be 
discussed latter on. We first summarize here some well 
known properties of the dimerized Heisenberg model (in 
the following we set J = 1 ) . 

Systems with 5 = 0. A spin gap smoothly opens at 
J 2c ~ 0.2412 (Ref. ^^HQ accompanied by a spon- 
taneous lattice dimerization. The ground state is then, 
twofold degenerate. The Majumdar-Ghosh (MG) pointt£l 
(J2 = 0.5 and S = 0) is believed to correctly describe 
the physics when J2 C < J 2 < 1/2. The elementary 
excitations are depicted as solitons (s) and antisolitons 
(sjlJ, topological defects between the two ground state 
patterns. No low energy ss bound states occur in this 
modelij and a continuum of magnetic excitations begins 
exactly at A 00 = A 01 = 2A so i where A so i is the minimal 
energy to create a soliton. 

Systems with 5^0. For J2 > J r 2 C , a non-zero 5 lifts the 
degeneracy between the two ground states and creates 
a linear potentialj-that confines solitons and antisolitons 
into bound statesO. A ladder of ss bound states occurs 
below the continuum which begins at 2A 01 r .L£. when a 
higher energy ss pair decays into two pairsB't2l (see Fig. 
|l|). The number of ss bound states has to be determined 
numerically panxi appears to increase as the dimerization 
is decreasedH'EEl. Turning on a small 5 with J2 < J2 C cre- 
ates a quite different situation. The explicit dimerization 
opens a spin gap A 01 cx <5 2 / 3 transforming the massless 
spinons of the Heisenberg chains into massive solitons. A 
mapping to the massive Thirring model gives evidences 
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for the existence of exactly one triplet and one singlet 
bound states with a universal ratio A 00 /A 01 = \fi (Ref. 
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FIG. 1. Singlet and triplet states in the linear potential. 
The continuum begins with a S=2 state (in finite systems). 



In this paper, the Raman spectrum of dimerized 
Heisenberg chains is numerically investigated by Exact 
Diagonalization (ED) and Density Matrix Renormaliza- 
tion Group (DMRG) calculations. In section ||, we dis- 
cuss the possibility of observing ss singlet bound states 
using Raman spectroscopy on pure systems. In section 
III, we consider the effects of non- magnetic impurity dop- 
ing on the Raman spectrum. Lattice r elax ation is studied 
in section [V . Results of sections ||, [in] and IV are ob- 
tained for parameters close to the MG point, useful to 
gain insight for a qualitative understanding. Finally, the 
relevance of these results to CuGe03 are quantitatively 
discussed in section [v| in relation with recent experimen- 
tal results. 



II. PURE COMPOUND 

Within the Loudan-Fleury theoryil, the Raman oper- 
ator is written asQcj: 



Hr = yjSj.Sj+i + -fSi.Si+2}. 



(2) 



Here, 7 is a microscopic parameter, and terms propor- 
tional to the Hamiltonian have been omitted. The Ra- 
man intensity Ir can then be obtained: 



I r (lu) = |(*o|^|*„)| 2 ^(c - E n + E ), 



(3) 



where |\& n ) is a complete set of eigenstates with energy 
E n (n = corresponds to the ground state) and 8 e is the 
lorentzian of width e. The Raman operator Hr depends 
on the value of 7. However, varying 7 only makes minute 
changes for the quantities of interest at least for the pa- 
rameter values used in this paper. Hence, we shall use 
the operator Hr = J2i &i-Si+i in the following. 

A typical Raman spectrum obtained by ED on a closed 
ring of 28 sites is shown in Fig. ||a for J 2 = 0.5, S = 
0.05. These parameters are chosen for two reasons: i) 
to minimize finite-size effects, ii) to study the spectral 
weight of the two (possibly three) ss bound states that 



numerically have been identified for these parameters^. 
As can be seen in Fig. ^a, the Raman intensity is zero up 
to A 00 , where the first bulk bound state (BBS) occurs. 
This one is by far the most pronounced feature in the 
spectrum. The second bound state also has a significant 
weight while the third one is very weak and difficult to 
distinguish from the continuum without a precise analysis 
of finite-size effects. 

The relative weights of the peaks are plotted in Fig. 
||b for lattice sizes from 16 to 28 sites. Note that the 
slightly larger dimerization 5 — 0.1 used here reduces 
finite size effects but the number of bound states is low- 
ered from three to two. At large lattice sizes, the weight 
of the two first peaks is constant and clearly non-zero 
signalling their bound state nature in contrast to the 
third peak which vanishes in the thermodynamic limit 
indicating that it is part of the continuum. An equiva- 
lent analysis for parameters closer to experimental values 
(J 2 ~ 0.36, 5 ~ 0.014) (Ref. [n]) is difficult due to large 
finite size effects. However, the physics should be simi- 
lar and it therefore appears reasonable to interpret the 
experimentally observed sharp peaks in pure systems as 
singlet ss bound states. 
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FIG. 2. (a) Raman intensity calculated for a L=28 site 
ring with J2 = 0.5, S — 0.05 and a broadening e = 0.1J. 
The positions of the three singlet bulk bound state (BBS) 
from Ref. [u] are indicated by arrows. The dashed line shows 
the onset of the continuum, (b) Relative weight of the three 
lowest energy peaks in the Raman spectrum for J 2 = 0.5 
and 8 = 0.1 as a function of the inverse lattice size. The 
presence of two BBS (o, □) below the continuum (0) can then 
be inferred. 



III. DOPED COMPOUND 

We now study the effect of introducing in-chain non- 
magnetic impurities. CuGeOs can be doped with low 
concentrations of Zn. The spin Zn substitutes the mag- 
netic Cu 2+ and in a first approximation the effect of such 
an impurity is simply to break the linear chains. Hence, 
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we can model the effect of Zn doping by studying chains 
of different lengths with open boundary conditions. 

First we qualitatively analyse the case of odd length 
chains. One end of the chain begins with a weak link 
— S). The chain end will now create a linear po- 
tential capable of binding a soliton. The ground state 
is made of a . so] i ton-bound to the chain end by this lin- 
ear potentia£3'E3E3a. Other higher lying bound states 
in the linear potential may also been identified and we 
shall refer to these states as edge bound states (EBS). 
These EBS exist up to an energy A 01 above the ground 
state, where it becomes energetically more favorable to 
create a free ss pair. As the invariance translation is bro- 
ken, a pseudo-continuum due to the ss dispersion begins 
there. The lowest lying state of the pseudo-continuum 
should consequently be made of the solitonic GS plus a 
ss pair and have a spin 3/2 (the impurity bound soliton 
and the ss pair do not interact in the thermodynamic 
limit). A two-particle continuum can be observed when 
two ss pairs can be created, i.e. for energies larger than 
2A 01 . 




ss-ss 
continuum 

ss quasi-continuum 
•- S=3/2 

S=l/2 bound 
states 

Ground state 

FIG. 3. The potential well created by an open boundary. 
The four soliton-impurity bound states shown here corre- 
spond to the parameters of Fig. 

@-H (Ja = 0-5, 8 = 0.05). 



We now turn to a more quantitative study. The 
number of EBS depends sensitively on the parameters 
of the Hamiltonian and previous parameters such as 
Ji = 0.5 and 5 — 0.05 should be optimal to observe 
EBS, with small finite size effects. Fig. ^ shows in- 
finite size extrapolations of the energies obtained with 
ED or a DMRG treatment for these parameter values. 
Four soliton-impurity bound states below the continuum 
can be clearly identified. As discussed before, the energy 
of the soliton-impurity bound state has a finite upper 
limit and the onset of the pseudo-continuum is indicated 
by a spin 3/2 state. We have checked that its energy 
with respect to the GS is, in the thermodynamic limit, 
the lowest energy of a ss bound state, i.e. the spin gap 
A 01 = (0.5200 ± 10~ 4 )J. More explicitly, one finds (see 
Fig. | for the notation) A™ 13 ~ 0.20J, A™ p ~ 0.35J 
0.47J. 



and A^ p 




• o 

0.008 



0.010 



FIG. 4. Lowest lying S z = | (o) and S z = § (X) exci- 
tations for the J2 = 0.5 and 5 = 0.05 Heisenberg chain as 
a function of the square inverse of the chain length L (up 
to L=25) obtained by ED. The energy reference is the GS 
for each chain length. Dashed lines represent a (j?, -rs) fit 
and the extrapolations to infinite sizes are indicated. DMRG 
extrapolations are also plotted (♦). 



The wave functions (Sf) obtained in a DMRG calcu- 
lation shown in Fig. |^ clearly support the previous inter- 
pretation. The four lowest states (a,b,c,d) show a soliton 
bound to the weak bond edge, and the soliton moves fur- 
ther away from the impurity when its energy increases. 
In contrast, in the lowest spin-| state (Fig. Us), a ss 
bound pair can be clearly identified in addition to the 
solitonic GS seen in Fig. ||(a). 
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FIG. 5. (Si ) as a function of the site i for the for lowest 
spin-i (a,b,c,d) and the lowest spin-| state (e) for J2 = 0.5, 
S = 0.05 on a 51-site chain using a DMRG algorithm 
(m=128). In (e), a magnification by a factor 6 has been re- 
alized right to the dashed line. The 7 first sites of the lowest 
bound state (a) (A) and a soliton-antisoliton pair on a 44 site 
lattice (•) are also shown for comparison. 
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As Hr is a spin operator, transitions between the 
S=l/2 ground state and higher lying S=l/2 bound states 
can occur and should be identifiable in the Raman spec- 
trum as well-defined peaks at energies A™ p . The Raman 
spectrum of a 25 site chain with open boundary condi- 
tions is shown in Fig. [}|a (plain line) together with the 
one obtained for a pure 24 site system (dashed line). The 
low energy structure is clearly due to the presence of the 
impurities and are signatures of EBS. The high energy 
part is the bulk contribution and is very similar to the 
one of the pure system (see section [nj). The main effect 
of introducing impurities is then to create additional sin- 
glet lines below the singlet-singlet gap A 00 in .the Raman 
spectrum as has been experimentally observedH. The first 
of the peaks has an energy A™ 15 which we interpret as 
the energy difference between the ground state and the 
first excited bound state of an odd length chain. 

The contribution of the EBS to the total spectrum as 
compared to the bulk contribution should scale like the 
impurity concentration i.e. 1/L. Hence, in order to cor- 
rectly identify these states in the relative spectrum the 
weights are normalized to the impurity concentration. 
The normalized relative weights of the edge bound states 
are shown in Fig. |6|b for L varying from 11 to 25 sites and 
J2 = 0.5, S = 0.05. In this case it is known that four edge 
bound states exist. Hence, three transition lines should 
be observable in the Raman spectrum. As shown in Fig. 
||b the first two transitions between edge bound states 
have a non-zero weight. The third transition has a very 
small weight and is almost undetectable. 

The presence of the impurity also leads to changes 
in the BBS resonances. Since the translation invari- 
ance is broken, both the S=0 and S=l ss dispersions 
are reflected at momenta q n = nir/L. Consequently, the 
BBS peaks are split on an energy range corresponding 
to the respective ss bandwidths although they keep al- 
most all their spectral weight. New small satellites should 
then be observed in Raman experiments forming a quasi- 
continuum. The possibility of flipping the spin of the 
impurity bound soliton may lead in rather small systems 
to a peak at the onset of the S=l ss quasi-continuum. 
The weight of this process becomes rapidly negligible for 
L > 15 sites (see the weight of the fourth peak in Fig. ^b). 
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FIG. 6. (a) Raman intensity calculated for a L=25 site 
open (plain line) and a L=24 site closed chain (dashed line) 
with J2 = 0.5, 5 = 0.05 and a broadening e = 0.05 J. The po- 
sitions of the first two singlet edge bound state (EBS) A™ 11 , 
A^ p are indicated by arrows, together with the three BBS 
of Fig. ^. The dashed line shows the onset of the bulk con- 
tinuum, (b) Relative weight of the four lowest energy peaks 
in the previous Raman spectrum normalized to the impurity 
concentration as a function of 1/L. The presence of two EBS 
(o, □) can then be inferred. A third one (A) has an almost 
zero although non- vanishing weight. The fourth peak (o) cor- 
responds to the first S=l/2 state of the quasi-continuum. 



For even length open chains, there are two possible sce- 
narios. If the chain begins and ends with a weak bond, 
the system is equivalent to the addition of two odd chains. 
Consequently there is a soliton-impurity bound state at 
each chain end and its edge excitations can be easily in- 
ferred from the odd system ones. If both the chain edge 
bonds are strong no edge bound states exists and the 
Raman spectrum should be similar to the bulk spectrum 
for pure systems. Note that this latter case with strong 
bonds at. the edges is favored when lattice relaxation is 
allowedEJ. 

As a consequence, in agreement with what was pro- 
posed in Ref. ^5], the previous analysis indicates the pres- 
ence of low energy peaks below the spin gap in inelas- 
tic neutron scattering spectra, due to soliton-impurity 
bound states. 



IV. ROLE OF THE LATTICE RELAXATION 

So far, the underlying lattice was assumed to be frozen 
in a dimerized configuration. However, recent work sug- 
gest that lattice relaxation can be important and lead to 
qualitatively new physics. Such effects are particularly 
crucial when inhomogeneities are introduced in the sys- 
tem by, for example, applying a magnetic fieldcSEZI or 
when doping with impuritiescJ. 

In order to go beyond the previous treatment, we shall 
assume here that the lattice distortion can adjust locally 
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to the magnetic modulation of the soliton wavefunction. 
It has been shownESa that, in this case, the soliton ac- 
quires both a spin and a lattice component. The de- 
scription of this magneto-elastic soliton can be realized 
in term of a simple Hamiltonian, 

h = Yyo- + <y£4 +1 + j 2 ^4 +2 ] + \ K Y,% 

i i 

+K X j2(-m, (4) 

i 

including an elastic constant K and where the local dis- 
tortions Si are determined self-consistently to minimize 
the total energy (using e.g. an ED method). The effective 
interchain coupling K± introduced here is, in fact, related 
to the physical elastic constant K± by the self-consistent 
relation K± = K±S where S — |(<$;)b| is the dimerization 
order parameter in the bulk (pure system) . Note that the 
lattice relaxation which we consider here is purely static 
(the distortion Si varies only in space) . A study of quan- 
tum lattice fluctuations is beyond the scope of this work 
(for more details on this issue see e.g. Ref. ^9|). We then 
implicitly assume that the time scale of the spin fluctu- 
ations is significantly shorter than the one of the lattice 
(adiabatic approximation). In this case, the ground state 
|\&o) contains a spin part and a distortion {Sf}. Since the 
Raman operator acts on the spin part only and the lattice 
relaxation time is much larger than the spin relaxation 
time scale, the eigenstates \^ n ) will have the same dis- 
tortion {<5°}. Then, once l^o) is obtained with the self- 
consistent method as explained above, the distortion is 
"frozen" for the computation of the spectra Ir(u>) using 
Eq. @. 

The interchain coupling K± isjnrucial in order to create 
a soliton-impurity confinements. The origin of this con- 
finement is clear; in the absence of K± the lattice exhibits 
a solitonic distortion pattern and the point where Si van- 
ishes coincides with the maximum of the spin wavefunc- 
tion (Sf) as seen in Fig. [7|a. Hence, in the region between 
the impurity and the soliton the lattice dimerization in 
the chain is then out of phase with the order parame- 
ter of the bulk, {Si) b = (— l) l S. A finite value of K± 
then leads to an elastic energy cost proportional to the 
impurity-soliton separation. The new equilibrium posi- 
tion is then obtained when the elastic "force" equilibrates 
the magnetic pressure that tends to delocalize the soliton 
away from the wall as shown in Fig. |^b. When K± — > oo 
while the lattice constant K is tuned to maintain a fixed 
value of the bulk dimerization S (K <~ K±/S = Kx), the 
topological lattice defect is pushed completely towards 
the impurity end. The lattice becomes then perfectly 
dimerized and one recovers the previous model. 
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FIG. 7. Results obtained on a L=21 site chain for the spin 
polarization (Sf) (o) and the lattice distortion Si (•) as a 
function of the site i (the impurity is located on site 0) for 
J 2 = 0.5. (a) K = 7.7, K ± ~ 0; '(b) K = 11.7, Kx = 0.02 
(Kx ^ 4). 



We shall focus here on the modification of the Raman 
spectrum when Kx becomes finite i.e. when the lattice 
can relax as described above. As shown in Figs. ||, the low 
energy EBS are robust features in the Raman spectrum. 
Note that we have fixed the bulk value of = S ~ 

0.05 by varying the parameters K and Kx simultaneously 
so that a direct comparison with Fig. ^ can be made. The 
finite size scaling of the weights in Fig. ^3 unambiguously 
establishes the existence of well-defined 5-peaks which 
can be interpreted as soliton bound states. 

For decreasing mean-field coupling Kx important 
changes in the spectrum take place with respect to the 
uniformly dimerized chain; (i) the characteristic energies 
of the EBS resonances decrease and (ii) a transfer of 
spectral weight occurs from the lowest energy pole to the 
next one. The three spectra shown in Figs. @a correspond 
to Kx = 0, 0.1 and 0.2 and to (realistic) anisotropy ratios 
Kx/K = 0, 0.21 and 0.34 respectively. The energy shift 
of the two low-energy peaks to lower energy is clearly vis- 
ible (together with the transfer of weight) when the lat- 
tice becomes more and more anisotropic. Indeed, for Kx 
decreasing from oo ( "rigidly" dimerized chain) to (iso- 
lated ID chain), at fixed 5 = 0.05, the energy of the first 
resonance decreases from 0.20 J to 0.13 J while the second 
peak moves from 0.35 J to 0.24J. As the interchain cou- 
pling is varied, a large transfer of weight is observed in 
the spectrum. Indeed, almost all the low energy spectral 
weight is carried by the first resonance for a large Kx 
(with a relative weight in the first peak 0.040 vs. 0.0067 
in the second one) whereas this weight is mainly located 
in the second peak for a vanishing interchain coupling 
(with a relative weight in the first peak 0.00065 vs. 0.071 
in the second one). In this latter case, only one peak, the 
second one, can be distinguished in the spectrum (see the 
dotted line in Fig. Sa). 
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FIG. 8. (a) Raman intensity calculated on a L=21 open 
chain with an adiabatic lattice for J2 = 0.5 and a broadening 
e — 0.05J. Only the low energy region below the BBS ss 
excitation is shown. The dotted, dashed and full lines corre- 
spond to K± = 0, 0.1 and 0.2 respectively. The values of K 
(= 7.7, 9.7 and 11.7 respectively) have been chosen to fix the 
bulk value |{<5i}_g| = S ~ 0.05. (b) Relative weight of the two 
lowest energy peaks in one of the previous Raman spectrum 
(K± = 0.2, K = 11.7) normalized to the impurity concentra- 
tion plotted vs 1/L. The dashed lines correspond to fits of 
the data using the form a + fr£^E-LziZii. 



V. RELEVANCE TO CUGEO3 

So far all calculations have been performed for param- 
eters close to the MG point in order to minimize finite 
size effects and clearly identify the bound states. We 
now discuss the implications of our results for more re- 
alistic parameter values relevant to CuGe03. Param- 
eters such as J ~ 160 K, J 2 ~ 0.36 and 8 ~ 0.014 
have been inferred by ED from high temperature sus- 
ceptibility-juaasurements and the zero temperature spin 
gap valueO'EllEj. Raman spectra obtained by ED would 
have very large size effects for these parameters, making 
an interpretation very difficult. Indeed, previous studies 
for J2 = or J2 = 0.24 ha^e seen no indication of soliton- 
antisoliton bound stateQEd. Hence, we turn to DMRG 
techniques, first omitting possible lattice relaxation ef- 
fects. Although the associated spectral weights cannot be 
obtained, the different gaps can be precisely evaluated us- 
ing even open chains: A 01 / J ~ 0.1802, A 00 /J ~ 0.2529 
and A 02 /J ~ 0.3606 (see Fig. |a). 
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FIG. 9. Variation of (a) the singlet-triplet A 01 (o), sin- 
glet-singlet A 00 (□) and singlet-quintuplet gap 2 A 01 (A); (b) 
the gap between the ground state and the first A™ 15 (o), sec- 
ond A 2 mp (A) and third A 01 (□) excited state for J 2 = 0.36 
and 5 = 0.014 as a function of the inverse square chain length 
using a DMRG algorithm for open chains. 



Using the same parameter values, it appears that three 
impurity-soliton bound states (Fig. |io|a,b,c) occur for 
odd open chains. These three bound states have a very 
extended but still localized wavefunction, signalled by an 
exponentially decaying part of (Sf) far from the chain 
end. For instance, the exponential part of the ground 
state (Fig. [Tc|a) is ~ cxp(— r/6.8). In contrast, the 
first S=3/2 state shown in Fig. nOp is made of a soli- 
ton close to the impurity (identical to the ground state 
(a) in the thermodynamic limit) plus a soliton-antisoliton 
pair which extends all over the rest of the chain. Note 
that the last bound state appears immediately below the 
first S=3/2 state. This is clear for short and intermediate 
lengths, however, this state could possibly become part 
of the ss quasi-continuum when L —> 00. Our data seem 
to indicate that this is not the case. 
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FIG. 10. DMRG results for (Sf) as a function of the 
site i for the three lowest S=l/2 (impurity-soliton bound 
states: a, b, c) and the lowest S=3/2 state (onset of the 
quasi-continuum: d) on a 51 site chain beginning with a weak 
link for J 2 = 0.36, 6 = 0.014. 
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The two gaps associated with these three impurity- 
soliton bound states have also been evaluated by DMRG 
(Fig. |b): Ai mp /J ~ 0.116, and A™ p /J ~ 0.175 
(just below the quasi-continuum) . The energy difference 
between the ground state and the third excited state, 
which has a spin 3/2 and is the first state of the quasi- 
continuum, gives approximately the spin gap A 01 (see 
Fig. ||b), as we have explicitly checked. Equivalently, us- 
ing J ~ 160 K~ 112 cm -1 , one obtains: A 01 ~ 20 cm" 1 , 
A 00 ~ 28 cm" 1 , of the same order as the experimental 
values ~ 17 cm" 1 and 30 cm" 1 respectively. The small 
deviation between theory and experiments may be at- 
tributed to other effects not included in this paper such 
as magnetic interchain coupling. For the doped systems 
our results are A" np ~ 13 cm" 1 , in surprisingly good 
agreement with the experimental result 15 cm" 1 (Ref. 
|8|). The ratio A™ 13 / A 01 is about 0.6, which is very sim- 
ilar to what has been found using a particle in a linear 
potentials. However, according to our analysis, another 
EBS should be observed at A™ 15 ~ 20 cm" 1 , contrary 
to the linear potential approach that predicts only onea. 
However, the spectral weight of this state might be very 
small. 

For an anisotropic lattice (see model (0)), another al- 
ternative explanation might have to be invoked. Indeed, 
in that case, significant lattice relaxation would shift the 
poles to smaller energies. One is then tempted to at- 
tribute the experimental resonance at 15 cm" 1 to the 
second theoretical EBS. In this scenario, another EBS is 
expected with a small weight at energies below 10 cm" 1 
A direct comparison with the experimental spectrumEl 
may corroborate this prediction but the signal noise pre- 
vents a definitive conclusion. Furthermore, detection of 
resonances in Raman scattering experiments at very low 
energy is difficult due to Brillouin scattering. 

In conclusion, we have shown that ss bound states 
should be easily identifiable in experimental Raman spec- 
tra on pure systems. Consequently, it is reasonable to 
interpret the experimentally observed line at 29.7 cm" 1 
as a ss bound state. Doping with chain breaking impu- 
rities should give rise to a number of edge bound states 
with non-negligible Raman weight. At least one such 
transition has already been experimentally observedu in 
CuGeOa. 
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time on the C94 and C98 CRAY supercomputers and 
rs6000 workstations. We also acknowledge funding from 
the ECOS-SECyT A97E05 programme. J.R. thanks 
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